Numerical generation of a vortex ring cascade in quantum turbulence 
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A pair of perturbed anti-parallel vortices is simulated using the three-dimensional Gross-Pitaevski 
equations. The simulations show decay of the vortex line length along with the depletion of local 
kinetic energy in a manner analogous to the decay of kinetic energy in the viscous, Navier-Stokes 
equations, even though the governing equations are Hamiltonian and energy conserving. The stages 
of development include vortex line stretching, conversion of kinetic energy into of the interaction 
component, the generation of vortex waves, multiple reconnect ions along the waves leading to the 
creation of vortex rings, kinetic energy transfer to high wavenumbers, evidence for energy transfer 
into phonons and finally the emission of the vortex rings. Some of this has been seen before, 
but not in a single calculation, which demonstrates the existence of a mechanism for an energy 
cascade that does not rely on statistical arguments. A four vortex example is given to demonstrate 
that these steps might be general. The wave generation and reconnection steps all depend upon 
interactions between distinct vortices, unlike vortex wave models where self-interactions along the 
vortices generate these effects. A variety of additional cases are noted whose analysis is in progress. 
So far, none have generated the dynamics that vortex wave models require to generate a cascade. 
Further cases are being pursued. 

PACS numbers: 47.37.+q,47.27.De,47.32.C-,67.25.dk 



I. BACKGROUND 

Despite the absence of viscosity, experiments have re- 
peatedly shown that superfluids exhibit resistance and 
depletion of turbulent kinetic energy in a manner simi- 
lar to the effects of turbulence in a classical fluid. [3T] 
made this quantitative, at relatively high temperatures, 
by finding that the vortex line density decays at a rate 
consistent with the kinetic energy decay of idealized clas- 
sical turbulence [11 . This result has since been confirmed 
at lower temperatures where the possibility of coupling 
with the more classical, normal fluid component is re- 
moved. These experiments include ^He for T < 0.2Tc 
^ , where Tc < 2.2mK, and more recent measurements 
in ^He at T < 0.5K [34]. Besides the similarities to decay 
in classical fluids, superfiuid spectral measurements yield 
classical k~^^^ power laws, both for velocity [23 and for 
a signature of the vortex line density [26 . 

Why would an inviscid, Hamiltonian system behave 
in this manner? Unfortunately, our understanding of 
classical turbulent decay is so poor that little insight 
is provided from that quarter. We do not even know 
if the equations we all believe underlie classical turbu- 
lence, the Navier-Stokes equations, are well-posed. What 
is observed in high resolution numerical simulations of 
the Navier-Stokes equations is that the stretching of fi- 
nite diameter vortices is crucial, there is a k~^^^ energy 
spectrum consistent with experiments and ultimately the 
kinetic energy is converted into heat by viscous dissipa- 
tion. Most theories explaining the cascade are statistical 
and are not based directly upon either the Navier-Stokes 
equations, or the vorticity equation. 

In this paper, a series of steps towards a quantum tur- 
bulent state with local kinetic energy depletion will be 
demonstrated in relatively modest three-dimensional nu- 
merical simulations of the Gross-Piteavskii equations, the 



nonlinear Schrodinger equation most relevant to quan- 
tum fluids. 
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Unlike in classical fluids, the velocity does not appear in 
these equations. Only through the Madelung transfor- 
mation does velocity appear and can defects in the wave 
function be interpreted as the equivalent of fluid vortices. 
Crucial differences with classical vortices are that quan- 
tum vortices are infinitesimally thin and all have the same 
value of circulation. Unlike the Navier-Stokes equations, 
these equations are a Hamiltonian system and conserve 
their global energy and are regular for all time. That is, 
no higher-order derivative of the wave function can blow- 
up (have a singularity) in a finite time. Despite these 
difference, reconnection of vortices occurs in the inviscid 
Gross-Piteavskii equations, which in classical systems is 
a purely viscous effect. 

The primary initial state to be discussed consists of 
two isolated, anti-parallel quantum vortices initially sep- 
arated by several core diameters (the healing length) with 
symmetric three-dimensional perturbations. It is found 
that these vortices transform themselves into a flow con- 
sistent with most of the observed properties of quantum 
turbulence. In particular, the kinetic energy spectra in- 
creases in slope to be nearly -5/3 and kinetic energy is 
removed from the interaction region, such that locally it 
appears as if the energy is decaying. It is found that wave 
dynamics and stretching on the vortex lines play comple- 
mentary roles in generating vortex reconnection. Neither 
by itself seems capable of generating the final state. 

The most important new dynamics to be identified wil 
be a mechanism by which an energy cascade to small 
scales can be generated from the interaction of only two 
vortices. The required steps are: 
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• Vortex stretching generates a vacuum between the 
points where the vortices are most strongly at- 
tracted. 

• There is an initial reconnection across this vacuum. 

• Curvature and torsion on the new vortex line cre- 
ated during the stretching process generates a sec- 
ond reconnection event and the release of a vortex 
ring. 

• During these phases, the vortex line length in- 
creases and kinetic energy is converted into inter- 
action energy.. 

• Waves form further out along the original vortex 
lines, then deepen. 

• There are further reconnections leading to a state 
dominated by vortex rings. 

• A k~^^^ kinetic energy spectrum forms at high 
wavenumbers in ky. 

• The interaction energy is released from the origi- 
nal interaction region, either by breaking down into 
phonons or by the rings propagating to the bound- 
aries. 

The paper will first present the equations, conservation 
properties, the numerical method and the initialization 
procedure. The discussion of the results will start by 
describing the evolution of the vortex topologies in phys- 
ical space. Next the development of the high wave num- 
ber spectra is discussed. Finally, the mechanism for the 
depletion of kinetic energy from the initial interaction re- 
gion will be presented. Several other cases are presented, 
including a four ring case that support many aspects of 
the scenario suggested by the anti-parallel case. 

A. Simulating superfluids 

To simulate a nearly zero-temperature superfluid, the 
following properties must be satisfied. 

• Because a quantum fiuid has no dissipation, its 
Hamiltonian, composed of gradient and interaction 
energies, should not decay in time except through 
additional physics or through interactions with the 
boundaries. 

• A quantum fluid is compressible, with a complex 
relation between the density and the pressure. 

• A quantum fluid is irrotational except along in- 
flnitely thin defects in the wave function known 
as quantum vortices and these all have the same 
quantum of circulation about their cores. 

The differences with classical vortices are: 



— In a classical fluid, vorticity is distributed uni- 
formly in space and the circulation about vor- 
tex cores depends on the initial condition. 

— For an ideal classical fluid, these values of the 
circulation are constant along Lagrangian tra- 
jectories. But in a viscous fluid the circulation 
changes in time. 

Simulations of the Gross-Piteavskii equations with 
good numerics and adequate resolution should satisfy 
these properties. However, these equations must be dis- 
cretized on a mesh and the numerical method must be 
remain accurate even if the density goes identically to 
zero in vortex cores. Until recently, the numerics seemed 
incapable of achieving these requirements. For this rea- 
son, vortex methods are often applied. However, vortex 
methods, whether Biot-Savart or the further reduced lo- 
cal induction approximation (LIA) have properties that 
might not be consistent with the Gross-Piteavskii equa- 
tions. Foremost are: 

• In LIA vortices do not stretch, so the vortex line 
length cannot grow. 

It will be demonstrated here, as have earlier Gross- 
Piteavskii calculations [2*, 'ST , that the line length 
grows just prior to reconnection, indicating that 
vortices can stretch. And subsequent to reconnec- 
tion, the line length decreases. 

• One expects, and these calculations conflrm, that 
reconnection can lead to waves, twisting and fur- 
ther reconnections. 

However, both LIA and Biot-Savart calculations as- 
sume sharp reconnections, which at least in the 
orthogonal case [18 , leads to exaggerated twist- 
ing and self-crossings that are not seen in Gross- 
Piteavskii calculations (TJ [15] . 

This will now be mentioned further to explain why 
these artifacts that do not appear in simulations of the 
Gross-Piteavskii equations nor in the experiments. 



B. Why not use a filament approximation? 

The flrst application of vortex methods to superfluids 
used the local induction approximation (LIA) [27]. In 
this approach, long-range Biot-Savart interactions are ig- 
nored, leaving only local dynamics along the vortex line, 
that is between neighboring points. Along a single line 
these equations can be transformed into a ID nonlinear 
Schrodinger equation [8 , which supports Kelvin waves 
and solitons. Both phenomena are observed in simula- 
tions of vortex lines, whether classical hydrodynamics, 
Biot-Savart fllament calculations, along GP vortices, as 
will be shown here, and perhaps most dramatically the 
twists that go up and down tornadoes. 
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While LIA explains some dynamical features for an 
isolated filament, dynamics due to interactions between 
vortices are excluded, including stretching and reconnec- 
tion. 



II. EQUATIONS, NUMERICS AND INITIAL 
CONDITION 



The advantage of considering only these two compo- 
nents of the energy is that this decomposition also applies 
in Fourier space and by examining the interaction terms 
in detail in Fourier space, one can more fully understand 
how a cascade might develop in that space. This is devel- 
oped further in the section on spectral analysis and will 
be the subject of a subsequent paper. 



The equations simulated are a form of the standard 
higher dimensional focusing nonlinear Schrodinger equa- 
tion for a complex wave function ip{x) wit h a cubic non- 
linearity: The Gross-Pitaevskii equations (I.l). This is a 



hard-sphere approximation to the full equation of state: 



(nV2m)VV(ic)- 



(II.l) 



Following [2j, the following choices are made: V{\x — 
x'\) = 0.b6{x — x')^ the chemical potential Ey = 0.5 and 
h and m are non-dimensionalized to be 1. While only 
an approximation for a superfluid, the hard-sphere ap- 
proximation is believed to be valid for a Bose-Einstein 
condensate and is believed to qualitatively incorporate 
the most important physics for quantum turbulence in a 
very low temperature superfluid. Therefore, while ideal- 
ized, it is physically relevant. 

The outer boundary condition for a BEG in a trap 
would have the wave function ip and density go to 
zero at the boundaries, that is Dirichlet. For a super- 
fluid in a container, it is better to use no-flux or Neu- 
mann boundary conditions, which is used here both as 
the outer boundary and as a way to generate symmetric 
initial conditions about the point of interest, as described 
below. 

Applying the usual Madelung transformation for quan- 
tum wave functions, if ip = ^^ex.^p{iO) then the velocity 
is 



These equations conserve mass: 



M 



{11.2) 



(11.3) 



and a Hamiltonian 



H=l J dV[\/^lj- VV^^ + 0.5(1 - (II.4) 

where ip^ is the complex conjugate of ip and 

K^^ = ^\V^f and Ei{x) = ^{l-\^ff (II.5) 

i^v^ is the kinetic or gradient energy and Ej is the in- 
teraction energy. 



A. Numerical method 

The numerics is a standard semi-implicit spectral algo- 
rithm where the nonlinear terms are calculated in physi- 
cal space, then transformed to Fourier space to calculate 
the linear terms. In Fourier space, the linear part of 
the complex equation is solved through integrating fac- 
tors with the Fourier transformed nonlinearity added as 
a 3rd-order Runge-Kutta explicit forcing. The domain is 
imposed by using no-stress cosine transforms in all three 
directions. 

This scheme has been used successfully for another 
ideal equation, the incompressible 3D Euler equations 
[6 where the timestep can be chosen using the Courant 
condition and the kinetic energy is conserved to within 
machine accuracy. Because the velocity has no role in the 
integration of the CP equations and because the nonlin- 
earity is cubic, not quadratic, it is not possible to choose 
the timestep in a similar way to conserve the total energy. 

Instead, a resolution dependent constant timestep is 
used that is chosen to be small enough that the results 
converge in a manner consistent with the Runge-Kutta 
algorithm. A measure of the accuracy of the method is 
to track the Hamiltonian over time, which grows slightly 
with time, roughly as {St)^. Over the course of the 
highest resolution calculation reported here, with 20,000 
time-steps, it grew by a factor of 0.01. 

The boundary conditions serve two purposes. Because 
the initial condition has two symmetries, it incorporates 
those symmetries into the numerics at the appropriate 
boundaries. Secondly, on the other boundaries the cosine 
transforms provide impenetrable, stress-free or free-slip 
boundaries, which are a reasonable approximation to the 
true boundary conditions in a superfluid. If more realistic 
conditions were required at the outer boundaries, extra 
terms could be added that could extract energy from the 
system. However, the results presented here demonstrate 
that even for the ideal equation there are effects that can 
mimic dissipation without relying upon any additional 
dissipative terms. 



B. Initial conditions 

The goal for initialization is to provide a smooth initial 
condition along a curved vortex. Due to the symmetries 
imposed, one could initialize with only one vortex. How- 
ever, it was soon realized that this resulted in disconti- 
nuities of higher-order derivatives of the wave function at 
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Time=0 £ =0.05 Density q: Y=0 Symmetry plane. Time=0 




FIG. 11. 1: Zooming in on the initial T — 0, vortex. Left: 
Three-dimension isosurface with p = .05. Right: Density 
contours through the y — symmetry plane of maximum 
perturbation. Recall that this is only one half of one of two 
mirrored vortices. A complete geometry for a similar initial 
condition appears in [6]. 



Time=1 p0.05 Y=0 Symmetry plane. Time=1 




FIG. II. 2: Early evolved time, T = 1. Left: Three-dimension 
isosurface with p = .05. Right: y = 0, p contours. First, as in 
a classical fluid, due to vortex stretching, the isosurface thins, 
most evident for the segment near y — z — 2. Near ^ = 
at this time, the stretching has pulled suflicient density away 
that an extended gap of p ~ is opening up. Reconnection 
will occur across this gap. 



the boundaries, and as soon as the cosine decomposition 
was imposed, waves appeared throughout the domain. 
These discontinuities were removed by using multiple im- 
ages from outside the domain of the initial condition and 
a Fourier smoothing. Earlier work [15] used three sec- 
ondary images. In this work, 23 secondary images are 
used. 

Another feature that is required for a smooth initial 
condition is that the density profile be both smooth and 
respects any expected analytic balance between centrifu- 
gal forces and pressure about vortices. The most sin- 
gular part of the initialization procedure is around the 
vortex cores, so to be on the safe side the initialization in 
the vortex cores was chosen smoother than might have 
been necessary. The final step for obtaining a wave- free, 
smooth initialization is to filter the highest wavenumber 
with a exp(— a/c)^) filter. Following experience gained 
initializing for the incompressible Euler equations [13 , 
the final filter should be used only to remove any final 
extraneous noise, and not be the primary means by which 
inconsistencies are removed. A measure of how successful 
all the other steps are at generating a wave- free, smooth 
state is how small the filter is needed, a = .002 when 
max(/c2) = 5 was found to be sufficient and in the worst 




FIG. II. 3: Leading up to reconnection, a kink appears at the 
transition between where the collapsing highly twisted vor- 
tex meets its straight continuation. Out of this kink, waves 
propagate outward along the vortex. After reconnection, the 
waves don't simply propagate. They deepen until a second 
reconnection occurs and a vortex ring separates for the origin 
vortices. 




FIG. II. 4: By foreshortening the y direction, the developing 
twisting can be better illustrated. 



test case reduced the pre-filter Hamiltonian by 2% and 
the density changed by .02%. 

The basic density profile about the singularities in the 
vortex cores obeyed an analytic profile of the form 

\^\ = Vp = r'^Nr^ + «o (11.6) 

While it is based upon a Fade approximate to an ideal 
2D quantum vortex [3 , which has m = 1 and ao = 2, for 
these calculations m = 2 and ao = 2 was chosen to ensure 
that the wave function around the vortex singularities 
was as smooth as possible initially. Once the calculation 
started, this profile quickly adjusted to m = 1 profile. 

The initial profile function was applied approximately 
perpendicular to the trajectory of the vortex lines, and 
not perpendicular to the 7/-axis as was done for an Euler 
calculation [6 . The reason for the change is that with the 
former procedure, there were large (> 50%) variations in 
the density above the background of p = 1 in the initial 
condition. Using the new method the maximum excess 
density peaks were reduced to p < 1.1. 

Two other differences with the Euler simulations are 
the trajectory of the vortex and the size of the domain 
in the vortical direction. Both of these modifica- 
tions were dictated by the formation of post-reconnection 
Kelvin waves on the vortices. For a smaller domain these 
would propagate to the ^/-domain wall and reflect back by 
the time only one vortex ring had separated off. To see 
the multiple rings reported required that the domain wall 
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be set back further, and that there be no perturbations 
to the initial anti-parallel trajectories beyond a short dis- 
tance in y. So rather than the nested cosine trajectory 
used for Euler [13 , the new trajectory path is: 



'cosh((^/(5^)i-8 



1,1,0 



(IL7) 



with 8x = —1.6 and Sy = 1.25 in di Lx x Ly x = 
Stt X Wtt X 47r domain on a 128 x 512 x 64 mesh. The 
power 1.8 on the normalized position y/Sy was chosen so 
as to localize the perturbation near the ^ = symmetry 
plane. This new trajectory assisted in minimizing the 
maximum excess density spots. 



III. RESULTS: STRETCHING, 
RECONNECTION, WAVES, RINGS 

The first three figures are designed show the stages in 
development before vortex rings begin to form. Figure 
|II.1| shows two early two isosurfaces of density for the 



initial condition chosen. Fig. II. 2 shows the initial vortex 



stretching. And Fig. II. 3 shows the evolution from this 
initial condition to a state where the vortex waves have 
started to refiect off the external wall. 

a. Stretching That vortices stretch in the Gross- 
Piteavski equations was first suggested by [24] using spec- 
tral properties. Using a box counting method [21] pro- 
vided evidence for both the growth and decay of vortex 
lines. How a pair can stretch is demonstrated here by 
Fig. |II.2| and quantified in Fig. |IV.1| For a superfiuid. 



vortex stretching changes the density between the vortex 
cores. First drawing denser fiuid in and making the re- 
gion of p ^ thinner, then as the two vortices approach 
one another, creating a vacuum. 

b. Vacuum and reconnection The creation of the 
pre-reconnection vacuum starts on the x — z symmetry 
plane then extends to the x — y dividing plane at t = 1. 
It is through this vacuum that reconnection occurs. The 
creation of the vacuum results in a transformation of ki- 
netic energy into interaction energy. 

c. Waves Following reconnection and the creation 
of new vortices with sharp curvature, Kelvin waves are 



generated along the new vortices, as shown in Fig. II. 3 



Similar waves were seen using the local induction ap- 
proximation following a wall-induced reconnection event 
[28] and have been explained using LIA [I8l [32]. The 
self-crossing of the vortex lines noted there bears some 
resemblance to the deepening of the vortex lines following 
reconnection reported here. 



A. Generation of rings 

The vortex waves continue to deepen until they recon- 
nect again across the x-y dividing plane in Fig. |II.3[ This 
is best illustrated by the T = 4.5 isosurface in Fig. II. 3 



Closer examination reveals a series of these kinks is de- 
veloping for y > 5 as the vortex waves become a series of 
zig-zags. A second reconnection at the kink near y = 5 
results in a ring separating from the paired filaments. 

After the second reconnnection and separation of the 
first vortex ring, a succession of reconnections quickly 
occur between the wavy vortices, that is between the il- 
lustrated vortex and its mirror across the x — y dividing 



plane, which lead multiple rings. Fig. |III A| shows the iso- 
surfaces and x — y dividing plane cross-sections of density 
to show where these subsequent reconnections occur. 

Note that each successive ring that forms has a smaller 
radius than the previous ring. Because the quantum 
circulation about the vortex cores of each vortex ring 
is identical, this implies that the propagation velocity 
V ~ T/R of each ring increases, and the separation 
between the rings increases. This implies they will not 
interfere with each other and can therefore freely leave 
the local system. 

Why are the successive rings each smaller than the 
previous one? Could it be an artifact of the particular 
initial condition of two nearly anti-parallel vortex lines? 
The four vortex case mentioned below demonstrates a 
state similar to this when the tangle is formed. Higher 
resolution and further analysis will be necessary to de- 
termine which secondary properties associated with the 
formation of multiple rings, exist in this calculation as 
well. 

Where, in physical space, could a cascade be form- 
ing that the spectral analysis reported next finds? Since 
the rings separate with time, it cannot be forming from 
subsequent ring-ring interactions and further reconnec- 
tion. Instead, the ongoing analysis suggests that each 
spectral cascade step is associated with the formation of 
each new ring. That is, each new and smaller ring pulls 
energy into its smaller length scale and corresponding 
higher wavenumber. 



B. Spectra 

Before spectral results are presented, some caveats are 
in order. Spectra are usually taken from large statistical 
ensembles, which allows them to smooth out transients. 
For time-developing calculations such as this calculation, 
smooth spectra are not expected. The most we can ex- 
pect is to see a progression with time of energy to smaller 
or larger wavenumbers. Therefore, the spectral slopes 
can be only approximate. Four spectral plots at succes- 
sive (but not equally spaced) times are shown. 

What might be significant is that the kinetic energy 
moves to higher wavenumbers in kx and , the directions 
perpendicular to the direction of the original vortex lines, 
and perpendicular to the direction of propagation of the 
final vortex rings. This would be consistent with the 
cascade developing through the formation of smaller and 
smaller vortex rings. 

The interaction energy spectra are equivalent to spec- 
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Tirae=12 c in dividing plane t=1 .2, case: vcosh3 




FIG. III.l: Isosurfaces [p < 0.05) and density contours in the 
X — y dividing plane between the two original vortices. Times 
are t = 12 and 48. Spectra at t = 12, 26 and 66 are shown 
below. 

tra of the density fluctuations, which have been measured 
in a superfluid [26 . The interaction energy spectra in 
these simulations are steeper at all times than -5/3, sug- 
gesting that Ej is being transferred to small wavenum- 
bers. 



IV. ENERGY DECAY 

A final observational property of quantum turbulence 
that this simulation can address is how dynamics analo- 
gous to viscous dissipation in classical turbulence could 
appear in an inviscid Hamiltonian system. What is ob- 
served experimentally is the decay of the vortex line 
length. Or more exactly the density of deficits in the 
vortex cores, which scatter the second sound waves or 
ion beams used as probes Assuming the lines are very 
thin, then the density of these deficits is proportional to 
the length of the vortex lines. 
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FIG. III. 2: Top frame: The initial one-dimensional 
K^^{kj) — spectra in {x,y,z) are shown, along 

with one interaction energy spectrum to demonstrate that 
initially the two components of energy are not strongly corre- 
lated. The initial spectra show resonances associated with the 
spacing, diameter and the perturbation chosen. This is espe- 
cially true for kz , where there is a bump at high wavenumbers 
that probably represents strong gradients within the vortex 
core which is oriented perpendicular to x and z. 
Bottom frame: By t — 12, just after the first reconnection, 
the K\/^ and Ei components follow roughly the same scal- 
ing in all three directions and to a rough approximation obey 
. A strong bump at high wavenumbers still persists from 
the initial condition for kz . While this looks like a bottleneck, 
it is only an artifact of the initial condition. 



A. Assumed relation between line length and 
kinetic energy 

This density of lines is usually interpreted as a sign of 
the decay of kinetic energy by the following argument: 

• Vortex line length in superfluid turbulence is ob- 
served to decay as /d. = t-3/^. 



• If each quantum vortex line segment is viewed as an 
element of vorticity cj, this implies that j dVuj'^ ^ 



. In classical turbulence, / dV.^ is known as the en- 

strophy, which (multiplied by viscosity) is propor- 
tional to the dissipation rate e of the kinetic energy 
E. 

So if dE/dt = e - t'^ , then E - 

• In classical, homogeneous isotropic turbulence in a 
periodic box the following is observed: 

- Energy: E(t) - t''^ 

— [TT]; but originally due to S. Patterson. 
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FIG. III. 3: Top frame: All spectra at t = 48 for reference. 
Bottom frame: By t — 48, for K\/^/J there is a long region 
in ky where the K^^ spectrum are the order of k~^^^. The 
interaction energy spectrum Ej is still dominated by a 
slope. Spectra in the other directions have similar trends but 
are less distinct. 



• Note that this decay law is never seen experi- 
mentally, as in all classical experiments there are 
boundary layers. It has only been seen in sim- 
ulations with either periodic or Neumann bound- 
ary conditions. As explained above, the Neumann 
boundary conditions used in the present calcula- 
tions are a reasonable first-order approximation to 
the true experimental superfluid boundary condi- 
tions. So these ideal classical simulations are rele- 
vant. 

• This is the basis for the claim that the observations 
of the decay of vortex line density ^[ |34] in super- 
fluids are related to the decay of kinetic energy in 
classical turbulence. 

When originally found in counter-flow at higher tem- 
peratures [31], the proposed explanation was transfer of 
energy to the normal fluid component through mutual 
friction. This would have been consistent with the ex- 
planation for relaxation in rotating bucket experiments 
going back to the 1950s. However, it has been pointed out 



that for this mechanism to work in a counter-flow experi- 
ments, there needs to be sufficient time for the superfluid 
and normal fluid vortices to equilibrate (Schwarz, private 
communication, about 1995). 

The more recent, lower temperature experiments [5l 
l3^ can completely rule out the normal-fluid coupling ex- 
planation. So how else could energy be removed from the 
superfluid? 



B. How could energy be depleted? 

All explanations center upon the boundaries. Unlike 
the ideal boundaries of these simulations, the physical 
boundaries in a superfluid cell are not ideal and will re- 
move energy, although not through viscous boundary lay- 
ers as in a classical fluid. Therefore, the question of the 
origin of superfluid decay reduces to finding a sufficiently 
efficient energy transport from the interior to the bound- 
aries to explain the observed decay of vortex line length. 

Three mechanisms have been proposed. 

1) Quantum vortex lines could reconnect, form vortex 
rings, which then propagate out of the superfluid 
tangle [7|. Some role for this mechanism is found 
here. 

2) Linear waves, or phonons could be generated and 
propagate to the boundaries if their production is 
strong enough. The importance of such a mecha- 
nism is demonstrated here, but not production of 
phonons by vortex motion. What is observed here 
would be consistent with the production of rarefac- 
tions and phonons in calculations of the collision of 
two vortex rings [20] and other configurations [2]. 

3) Recently, there have been several proposals based 
upon the propagation of waves on the vortices, 
which Fig. |II.3| shows can develop on Gross- 

These wave theories could pro- 



Pit eavskii vortices, 
vide the decay mechanism if they transport en- 
ergy sufficiently rapidly to their ends, presumably 
on the boundaries, or if there is a wave cascade 
to small scales along the filaments, at which point 
the smallest waves could generate phonons, which 
would then move the energy to the boundaries. 

For the Kelvin wave explanation to work, the cas- 
cade rate must be sufficiently efficient. Traditional 
wave cascade models have not been able to achieve 
this efficiency, as pointed out by [2?. Alternatives 
have been proposed [16, J/7 based upon adding ad- 
ditional physics and assumptions that would be dif- 
ficult to verify experimentally. 



C. Kinetic energy depletion in the calculation 

How can direct simulations of the Gross-Piteavskii 
equations help? The question is how to make the phe- 
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FIG. IV. 1: Estimates of the line length compared to changes 
in the interaction and kinetic energies. Followed by distribu- 
tions of the interaction energy with respect to the density at 
three times noted in the upper half atT = 0.5, T = 6 and 
T = 30. Analysis was done over the full domain (upper frame) 
as well as within quadrants defined in the ^/-direction, with 
only the first ^-quadrant plotted, which is where the interac- 
tions began. There is strong vortex line growth for T = 0.5—6, 
followed by a period for T = 6 — 25 of strong decrease in L. 
Distributions of the Ei with respect to density are shown at 
T = 0.5, T = 6 and T = 30 changes in the distribution of 
Ei during these stages. For T > 30 the global kinetic energy 
K\/^ grows again. This is associated with the accumulation 
of energy for large y and the rings colliding with their images 
at the opposite wall. In a physical cell, this energy would 
be absorbed by the outer wall. To show that kinetic energy 
and line length continue to decrease in the original interaction 
region, the middle frame shows K\/^ and L only in the first 
quadrant. The final dashed curve is a quantum equivalent of 
the Taylor microscale At- 



nomenological relationship of the line length decay to the 
kinetic energy decay more robust. With a calculation, all 
the components of the energy can be calculated directly, 
which we want to related to a consistent, robust algo- 
rithm of determining the vortex line length. 

Two methods have been proposed for determining the 
line length in calculations. Either by counting boxes with 
a property associated with the vortex cores [2T1. Or by 
making the assumption that the spectral integral e = 
'^k'^K\7^{k) = '^k^\i^{k)\'^ is proportional to the line 
length squared [24] by the arguments above under Sec. 

The analysis below shows that with the box-counting 



proposed here, the development of the anti-parallel calcu- 
lation might provide the first numerical evidence of when 
both line length and kinetic energy in the original inter- 
action region decrease. Elements of this can be found in 
earlier work, but either the box counting analysis is not 
compared to the kinetic energy [U |2T] or the spectral line 
length analysis is unphysical. That is e grows while ki- 
netic energy decreases [2^ , or in the present calculations 
e does not decay. 

The box-counting criterion chosen here is to a pick 
low density threshold Pc that will estimate of the volume 
occupied by the vortex cores. By div iding this volume 
by the area of the analytic profile (II.6) with m = 1, and 



choosing different density thresholds, consistency checks 
can be made. 

There are two primary stages to the time dependence 
of the line length. After a brief adjustment from the 
initial condition, completed by T = 0.5, the line length 
grows, which would be consistent with vortex stretching. 
The stretching stage lasts beyond the first reconnection, 
near T = 1.25, and on through the second reconnection 
near T = 4.5 when the first and largest vortex ring sep- 



arates off, as illustrated by Fig. II. 3 

Thereafter the line length decays. To understand the 
different stages, below the main figure there are three 
subplots showing distributions of the interaction energy 
Ej (II. 5) with respect to density at different times. If 



vortex stretching is occurring, one might expect there to 
be a strong growth in regions with very small densities. 
If the energy is in waves, then we might expect the distri- 
bution of Ej to be concentrated around p = 1. Around, 
not at, because for p = 1, = 0. Distributions of the 
density are not useful since their maximum is peaked near 
p = 1 for all times. The distribution at T = 0.5 is meant 
to demonstrate that initially, Ej is maximum near p = 1. 

During the period of vortex stretching there is a dra- 
matic growth in Ej, with most of the growth at small 
values of p, as demonstrated by the distribution for 
T = 6 when the line length L is maximum. Because 
Ei{x) = (1 — p)^, small p at a given x automatically im- 
plies a small value for Ej at a?, so what the distribution 
at T = 6 is telling us is that there is a large growth in 
the number of points with p « 0. Note that the strong 
increase in Ei during each of these stages is compensated 
for by a strong decrease in the kinetic energy 

Immediately after T = 6, L begins to decrease dra- 
matically while kinetic energy i^v^ continues to decay, 
which is compensated for by a continuing increase in the 
interaction energy Ej. At the end of this stage, there is 
a growth in large values of Ei on either side of p = 1, 
shown by the distribution at T = 30. This would be 
consistent the development of waves and the visualiza- 
tions of waves being emitted from two colliding vortices 
[U |20] . The double humped distribution relaxes gradu- 
ally for T > 30. 

A decrease in the global kinetic energy does not per- 
sist. Eventually interaction energy converts back into ki- 
netic energy. This could be associated with the structures 
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(rings) formed by the multiple reconnections reaching 
the outer ^/-boundary, or due to the existence of a wave- 
dominated state. Similar oscillations where observed in 
GP calculations with a symmetric Taylor- Green type of 
initial condition [24 . This would not happen in a real 
experimental device because the boundary conditions are 
not ideal. 

To mimic what decay in a real flow might look like, 
especially experiments that generate tangles far from 
boundaries [34^, the final two curves in Fig. IV. 1 show 
the growth, then decay, of the line length in the first y- 
quadrant (0 < < 47r), and a rescaled kinetic energy 
decay curve. When viewed in this way, the local region 
exhibits kinetic energy depletion due to transport out of 
that region in physical space after the global kinetic en- 
ergy begins to reappear. 




Orthogonal case vyz24r 
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V. ANALYSIS AND CALCULATIONS IN 
PROGRESS 

For this calculation to be accepted as a paradigm for 
the dynamics of kinetic energy decay in quantum tur- 
bulence, it must be demonstrated that the events seen 
here are not unique to this particular initial condition. 
The four ring case noted below would support this claim. 
What the anti-parallel case provides that more general 
cases cannot provide is a tool for addressing deeper ques- 
tions about the cascade dynamics, and in particular the 
relationship with the creation and destruction of physical 
space structures. What should be addressed? 



Nonlinear spectra: Formulation and cascade trends 

An unfinished piece of analysis is the Fourier spectra 
of the energy transfer. The equations have been formu- 
lated, and one property is that the kinetic energy and in- 
teraction energy must mediate the cascade of the other. 
There cannot be a cascade of just kinetic energy or inter- 
action energy. They must be coupled. Which directions 
do the two cascades go? The evidence so far suggests that 
there is a cascade of kinetic energy to small scales, while 
the conversion to interaction energy is through a cou- 
pling with much smaller wavenumbers (i.e. large length 
scales). 



More general initial conditions 

Upon presentation of these results at several meetings, 
one consistent concern was whether the dynamics ob- 
served here are an artifact of the particular strongly anti- 
parallel initial condition, or might be more general. 

One of the goals of these additional cases will be to 
test the validity of assumptions made by recent vortex 
wave models for the quantum turbulence cascade. So 
far, all wave growth and all reconnections seen are driven 



FIG. V.l: Two orthogonal vortex liness, the time dependence 
of the kinetic and interaction energies, two views just after 
reconnect ion. 



only by interactions between vortices that were originally 
distinct. There has been no evidence for any signficant 
wave growth or appearance of reconnections due to in- 
dividual vortex dynamics. Rather, the generation of Ej 
due to stretching and the phonons that are generated in 
the process seem to suppress all auto- vortex wave and re- 
connection processes. In particular there is no evidence 
for cascades of waves appearing on individual vortices as 
one model with higher order Biot-Savart terms predicts 
|T7]. And reconnections are never sufficiently sharp to 
generate the torsion needed to drive self-reconnections 
on vortices as proposed by another model ^16j . 



A. Orthogonal 

The case most frequently asked about is initially or- 
thogonal vortices. Test calculations were done and are 
essentially no different than what was first presented in 
1993 [15 . The reconnected vortices are twisted, but not 
enough to generate waves or any further self-crossings 
and reconnections. This is distinctly different from the 
LIA filament analysis of p!8^ , from which at least one ad- 
ditional one additional reconnection, would be predicted. 

This orthogonal case also resembles Case A from [4 for 
the experimentally visualized quantum vortices recently 
highlighted in in Physics Today j3Q] . 

Based upon simulations of the reconnection of classi- 
cal orthogonal vortices, one would have expected there 
to be self-reconnections. Further tests with the original 
vortices are intermediate angles need to be explored to 
determine if this result is general. 
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FIG. V.2: Two vortex rings combine via reconnect ion, then 
separate. 



FIG. V.3: Four colliding rings, the entangled state generated, 
the final relaxed state, and the time dependence of the kinetic 
and interaction energies plus a measure of line length in the 
inner region that contained the original vortices. In this case 
the measure of line length, the volume where p < 0.1, tracks 
the interaction energy Ei more closely than the kinetic energy. 



B. Two rings 



C. Four rings 



In the search for a structure with collapse properties 
similar to the primary anti-parallel case here, a number of 
colliding ring configurations were explored. One batch of 
tests used approximately anti-parallel pairs, all with their 
axes in the x — y plane, each with slightly different radii, 
different ellipticities and displaced horizontally from one 
another. 

All except one followed roughly the same scenario. As 
they approached, first there would be two reconnections. 
This would then form a roughly figure-eight configuration 
of the vortex lines. Then with further reconnections this 
would break down into two outgoing vortex rings. 

There are similarities to the dynamics of two collid- 
ing rings with the same radii, and displaced horizontally, 
reported by [20]. They do not get a figure-eight, but a 
more closed intermediate structure. However, they do get 
two rings coming out of the intermediate structure. The 
purpose of their paper was to demonstrate that phonons 
(or CP Kelvin waves) appear from the reconnections and 
possible evaporation of the smaller outgoing rings. This 
would be consistent with the anti-parallel case here. Fur- 
ther analysis of higher resolution calculations would be 
necessary to determine if this is happening for these test 
initial conditions. 

In one distinct case here, the combination of ellipticity 
and displacement of the vortices was sufficient that they 
reconnected at opposite ends of their ellipses, at the point 
where originally the vorticity was parallel. But first, one 
of the rings the vorticity twisted itself around so that 
the reconnection was still between approximately anti- 
parallel vortices. 



The best generalization of the anti-parallel calculations 
is to let four vortex rings collide, shown in Fig. |V C 



These temporarily combine to produce of tangle of wavy 
vortices, which breaks up into a variety of small rings 
that appear to evaporate (perhaps turning into phonons) 
and a large vortex ring. The large ring probably removes 
most of the original Hamiltonian. 

Three figures show the initial condition, the state of 
maximum entanglement, and the final state of a single 
large vortex ring. If the initial state is envisaged as two 
anti-parallel pairs, to get the entangled state it was nec- 
essary that the separation between the pairs be the order 
of their radii, but not overlapping. 

The fourth figure shows the time development of the 
kinetic and interaction energies. Again there is a strong 
initial transfer of energy from the kinetic to interaction 
components. The entangle time shown is roughly at the 
end of this process. After that, as structures start in- 
teracting with the boundaries, the interaction energy is 
converted back into kinetic energy. Line length analysis 
is needed as well as local analysis around the region of 
the initial interactions to determine if some of the trends 
observed in the anti-parallel case are repeated. 

Before this is done, this case needs to be rerun, first 
with much higher resolution to see the details of the en- 
tangled state, then with a combination of higher resolu- 
tion and larger domain to verify that the final state is 
not too strongly infiuenced by the boundary conditions. 

There are similarities to the four ring case investigated 
by [21 . The four rings in that case are all oriented per- 
pendicular to a common axis and have more symmetry 
than the case here. Resolution is also less than half that 
used here. Nonetheless, two primary wavy, vortices de- 
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velop, plus smaller entities that appear to break down 
into waves. 



VI. OUTLOOK 

For many years, simulations of first Navier- Stokes tur- 
bulence pj] have been able to generate k~^^^ spectra and 
turbulent structures, now getting 2 decades using up to 
4096^ mesh points [10 . More recently it has become pos- 
sible to do the same with the Gross-Piteavskii equations 

m- 

However, there are no clear relations between the spec- 
tra and structures in any of these cases. 

What would be particularly useful would be a single 
case, no matter how special, that could demonstrate a 
connection between vortex structures and the origin of 
the energy cascade. 

There have been special cases in the past showing a 
possible connection between structures and spectra, for 
example the Lundgren vortex [19 , but that case requires 
the existence of an imposed large scale stretching. More 
interesting would be a case that generates a flow of energy 
across scales that starts with individual structures. 

The case presented here provides such an example, and 
with time we should be able to demonstrate most of the 
underlying steps in more general cases similar to the fl- 
nal example given. Even then, the value of the special 
anti-parallel case is that it will be easier to identify the 
dynamics. 

Some of the dynamics discussed here have been noticed 
before as referenced above. This includes initial develop- 
ment where the length of the vortex lines grows, that 
is the existence of vortex stretching, the generation of 
waves on vortices, reconnection between lines and rings, 
and the determination of spectra. 

New relations between some of the known properties 
are noted here. For example, waves following artiflcial 
reconnections in the Biot-Savart and local induction ap- 
proximations have been seen, but waves on individual 
quantum vortices following reconnection have not been 
seen before. And from that, one can see is how those 
waves can lead to further reconnections, the formation 
of multiple rings, and an energy cascade. That is, this 
calculation provides, within the context of the equations 
for an ideal quantum fluid, a means by which waves on 
vortices can generate an energy cascade, a long-standing 
proposition for being the underlying source of the ob- 
served decay of quantum turbulence. And not requiring 
fllament approximations of any type. 

The results here suggest that even at higher temper- 
atures the decay of the quantum fluid occurs indepen- 
dently of the normal fluid. To determine this numerically 
would require a model that could couple the quantum 
and normal fluid components of a superfluid. Currently 
a good microscopic theory of the normal fluid does not 
exist, so phenomenological models would have to suflice. 
We would want a model where the coupling went to zero 



as T ^ [9]. 

The flnal new analysis shows how the transformation of 
kinetic energy into interaction energy is associated with: 
First, vortex stretching. Then the formation of rings, fol- 
lowed by a signature of waves appearing in distributions. 
The formation of waves does not seem to be due to forc- 
ing of the background quantum fluid by the vortices and 
their reconnections, but instead originates from interac- 
tions between the stretched, reconnecting vortices. 

VII. SUMMARY 

The analysis of the dynamics or reconnecting anti- 
parallel quantum vortices has identifled the following 
steps: 

• First, before vortex cores interact directly, Biot- 
Savart vortex dynamics of stretching, curvature 
and torsion dominates in a manner consistent with 
fllament calculations and simulations of the classi- 
cal, ideal Euler equations. 

During this phase, kinetic energy is converted into 
interaction energy by the vortex stretching and 
growth in the line length. 

• Once the cores begin to interact, the dynamics of 
the Gross-Piteavskii equations takes over, with re- 
connection developing in the vacuum that forms 
between the pair. 

• Following the flrst reconnection event, vortex waves 
are emitted with properties similar to waves in the 
local induction approximation. 

• These waves propagate down the initial vortex and 
deepen. 

• After they have deepened enough, secondary recon- 
nections occur and vortex rings form. 

• Near the time of the second reconnection, interac- 
tion energy pulled from the original kinetic energy 
is concentrated in the p ~ regions and spectra 
have a regime. 

• As the vortex rings fully separate, interaction en- 
ergy continues to grow, now it appearing more on 
either side of p = 1, indicating either waves or rar- 
efactions in the centers of the newly created rings. 

• During this stage, the high wavenumber spectra 
grow until spectra in two directions develop nearly 
-5/3 subranges. 

• All of these steps occur without any reduction in 
the global Hamiltonian. 

The parting question posed by these calculations is to 
flnd out how much of this could be carried over to dy- 
namics in the classical equations. One point of view could 
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be that despite the differences between the two systems, 
because the vortices in GP are weh-defined and there are 
no singularities, this might be a simpler system within 
which to identify some of the dynamics behind classical 
turbulence. To answer this, new classical calculations 
have begun with new initial conditions for the trajecto- 
ries of vortex lines based on the calculations here. 
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